A PDE for the joint distributions of 
the Airy Process 
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Abstract 

In this paper, we answer a question posed by Kurt Johansson, 
to find a PDE for the joint distribution of the Airy Process. The 
latter is a continuous stationary process, describing the motion of the 
outermost particle of the Dyson Brownian motion, when the number 
of particles get large, with space and time appropriately rescaled. The 
question reduces to an asymptotic analysis on the equation governing 
the joint probability of the eigenvalues of coupled Gaussian Hermitian 
matrices. 

The differential equations lead to the asymptotic behavior of the 
joint distribution and the correlation for the Airy process at different 
times t\ and ti-, when ti — t\ — > oo. 
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1 Main result 

The Dyson Brownian motion (see [2]) 

(Ai (£),..., A n (*)) el", 
with transition density p(t, /i, A) satisfies the diffusion equation 

dt 2^d\i d\i $(A) 



^ 1 3A~ 3A~ 



p 



with 



$(A) = A 2 (A) JJe- A '/ fl . 
i 

It corresponds to the motion of the eigenvalues (Ai(i), . . . , A n (t)) of an 
Hermitian matrix B, evolving according to the Ornstein-Uhlenbeck process 

dP „ , d 2 Id 



ftp /I f) z ] f) \ 

m^^+^m^^m^)' (L1) 

1,3=1 v J 

with transition density (c := e - */ 0-2 ) 



P(t,B,B) = Z- \, „, w2/2 e"^) 



Tr(B-cB) 2 

c 2jn 2 /2 



The -Bi/s in (1.1) denote the n 2 free parameters in the Hermitian matrix B. 
In the limit t — > oo, this distribution tends to the stationary distribution 

Z- 1 e~?' ItB2 dB = Z- 1 A 2 {\)]le-J d\. 



With this invariant measure as initial condition, one finds for the joint dis- 
tribution: 



P(B(0) 6 dB u B(t) 6 dB 2 ) = z-^^e-J^^W-* **». 

(1.2) 

Setting a = 1, the Airy process is defined by an appropriate rescaling 
of the largest eigenvalue A n in the Dyson diffusion, 



A(t) = lim V2n 1/6 (\ n { n - 1/3 t) - y/2n) , (1.3) 

in the sense of convergence of distributions for a finite number of t's. This 
process was introduced by Prahofer and Spohn [S] in the context of poly- 
nuclear growth models and further investigated by Johansson [Sj. Prahofer 
and Spohn showed the Airy process is a stationary process with continuous 
sample paths; thus the probability P(A(t) < u) is independent of t, and is 
given by the Tracy- Widom distribution 

P(A(t) < u) = F 2 (u) := exp (- J (a - u)q 2 (a)dc?j , (1.4) 

with q(a) a solution of the Painleve II equation, 

-fcJ 

- 9 e /= 171 for a /" oo 
g" = ag + 2g 3 with g(a) { W (1.5) 



s/—a/2 for a \ — oo. 

At MSRI (sept 02), Kurt Johansson posed the question, whether a PDE 
can be found for the joint probability of the Airy process; see [3]. The present 
paper answers this question, which enables us to derive the asymptotics of 
the large time correlation of the Airy Process. We thank Kurt Johansson for 
introducing us to this process. 

Theorem 1.1 Given t\ < t 2 , the joint 'probability for the Airy process 



V + x y — x, , .. . y + x ., . y — x 

h{t2-h) U —- ,2— ) :=logP L4(ti) < ?-—,A(t 2 ) < y —— 
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satisfies a non-linear PDE 1 in x, y and t — t 2 — t\, 



2t 



d*h 



dtdxdy \ dx 



d 



d\ (d 2 h d 2 h 



f d 2 h d 2 h 



dy)\dx 2 dy 2 ) \ dxdy ' dy 2 J y ' 

with initial condition 

y+x y-x\ ( y + x y-x . 

- logF 2 mm(— — , — — ) . 



;i.6) 



lim h I t: 

t\o V ' 2 



Conjecture For any fixed t > 0, i6l, the conditional probability satisfies: 
lim P(A(t) >x + z\ A(0) < -z) = 0. (1.7) 



Accepting this conjecture, we prove: 

Theorem 1.2 For large t, the joint probability admits the asymptotic series 

4/ X X rw NTW N F'(v) $ ( «, U ) + $ ( V , u) 

P(A(0) < < V ) = F 2 (u)F 2 (v)+ 21 ^ 2l ' + _L_lX_Ij ^ +0 I _ j , 



li 



wit/i «s the function (1.5)) 



$(u,u) := F 2 (w)F 2 (t;) 



oo \ 2 / r°° \ 2 
. q 2 da ] I / g 2 cfo: ] 

+ g 2 ( M ) U q 2 {v)-\(r q 2 day 



POO POO 

- da(2(v- a)q 2 + q' 2 - g 4 ) / g 2 cfa 
Moreover, the covariance for large t behaves as 



V 



where 



E(A(t)A(0)) - E(A(t))E(A(0)) = ! + £ + ... , 



c := 2 / / §(u,v)du dv. 



;i.9) 



Mn terms of the Wronskian {/(y), : = f'{v)g{y) - f{v)9'{v)- 
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Remark: The equation (1.6) for the probability 



/i(t 2 -*i;u,v) := logP(A(ti) < u, A(t 2 ) < t>), t = t 2 -h 

takes on the alternative form in the variables u and v, 

d t d 2 d 2 \ , <9 3 /i / d 2 h 8 2 h 8 2 h 

2 1 

2 dv V dv 2 dudv du 2 

d 3 h /&h d 2 h d 2 h 



dt^du 2 dv 2 ) du 2 dv^ dv 2 ~^ dudv du 2 



dv 2 du V du 2 dudv dv 2 
/<9 3 /i d d 3 h d \ / d ^ d \^ (1 10) 

\du 3 dv dv 3 du/\du dv) 



with initial condition 



limh(t;u, v) = logF 2 ( m i n (w, v )). 

This equation enjoys an obvious u v duality. 

The proof of this theorem is based on a PDE, which was obtained in pQ 
for the spectrum of coupled random matrices. 

In Spohn and Prahofer pose the question about the asymptotics of the 
covariant functions of A(t) and A(0) for large t. Moreover, in a very recent 
paper, Tracy and Widom [7] expressed the joint distribution, for several 
times t\, . . . , t m , as the exponential of a certain integral; its integrand involves 
traces of matrices, which satisfy a coupled system of non-linear ODE's. The 
quantities involved are entirely different and their methods are functional- 
theoretical; it remains unclear what the connection is between the two results. 



2 The spectrum of coupled random matrices 

Consider a product ensemble (Mi, M 2 ) G H 2 := H n x H n oinxn Hermitian 
matrices, equipped with a Gaussian probability measure, 

C n dM 1 dM 2 e -iT r( A^+M|-2 C M 1 A/ 2 ) ) (21) 

where dM 1 dM 2 is Haar measure on the product TC^, with each dM iy 

n n 

dMx = A 2 n {x) Yl dxid\J x and dM 2 = A 2 n {y) JJ d Vi dU 2 (2.2) 
i i 
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decomposed into radial and angular parts. In jl] , we define differential oper- 
ators Ak, Bk of "weight" k, which form a closed Lie algebra, in terms of the 
coupling constant c, appearing in (|2.ip . and the boundary of the set 

E = E x x E 2 := U[ =1 [a 24 _ 1; a 2i ] x U s i= ifoi-i, b 2i ] C M 2 . (2.3) 

Here we only need the first few ones: 

I i r d s d i \ f r d s d 



da-i dc 4— ' 9c 



(2.4) 



In PP, we prove the following theorem: 

Theorem 2.1 Given the joint distribution 

P n (E) := P(all(Mi- eigenvalues) G £i, all(M 2 - eigenvalues) 6 £2), (2-5) 

the function F n (c; a 1; . . . , a 2r , 61, . . . , 6 2s ) : = log -P n (-£-)? satisfies the non-linear 
third-order partial differential equation 2 : 

TIC I I TIC I 

S 2 ^!F n , BiA^n + -^—r } ~ < A 2 B 1 F n , A x E x F n + — — \ = 0. 



_2 _ 1 I 1 

ci J.4i Is ci JBi 

(2.6) 

Corollary 2.2 For E — Ei x E 2 : = (—00, a] x (—00,6], z't z's convenient to 
use the new variables x := —a + c6, y := — ac + 6. In these variables, the 
equation \2. 6}) for 

f n (c;x,y) :=\ogP n {E)=F n (q |^|) (2.7) 

ta£;es on the following form: 

d_ ( (c 2 - l) 2 g^ + 2ncx - n(l + c 2 )i/ 

dx \ (^ 2 -i)£fe + ^ 

«9 /(c 2 - l) 2 g^ + 2nc?/ - n(l + c 2 )x 

^ (^ 2 -i)Sfe+^ 



(2.8) 



2 in terms of the Wronskian {/, g}x = ^/-<? — f-Xg, with regard to a first order differ- 
ential operator X. 
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Proof: It is an immediate consequence of Theorem 2.1, upon observing the 
simple form of the differential operators A\ = d/dx and B\ = d/dy, when 
expressed in terms of x and y. ■ 



3 Proof of Theorem 1.1 

Taking into account the scaling in the definition of the Airy process (11.3(1 . 
and using the Ornstein-Uhlenbeck transition probability (jl.2j) . we compute 
the probability (setting c = e~ n lfit , with t = t 2 — ti) 



P (v / 2n 1/6 (A n (n- 1/3 t 1 ) - V2n) < u, V2n 1/6 (X n ( n - 1/3 t 2 ) - v 7 ^) < v) 



all ^-eigenvalues < -^(2n 1 / 2 + n 1 ^ 6 u) 



all £? 2 -eigenvalues < 4 5 (2n 1 ' /2 + 



x dP>idP>2 ■± I Tr(Bl+Bl-2cB 1 B 2 ) 

(l_ c 2)nV2 e 



all Mi-eigenvalues < ^j""^ 6 " 
all M 2 -eigenvalues < 



2n 1 / 2 +n~ 1 /S 

Z'- 1 dM 1 dM 2 e-^ M ' +M ^ 2cM ^\ 

(3.1) 



using the change of variables 

Bi 



Mi 



V0- - c 2 )/2 



The integral (|3.1(1 turns out to coincide with the statistics (|2.5() of the coupled 
matrix model with the Gaussian distribution (|2.1j) . Therefore, we set 



2U 1 / 2 + n-V« u 2^ + n- 1 / 6 



and 6 = A ; (3.2) 



Vl^ ^ 



c 



in formula (j2.8jl . via a; and Then, setting = n 1//6 , we now express x and 
y in terms of u, v and c = e~ n 1/a * = e~ t//fc2 , using the change of variables in 
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Corollary 2.2, 



x(c(t),u,v) = -a + bc = — l[2kr + T - [2k* + T c 



VI -c 



2/(c(t),tt,v) = -ac + 6 = ~^=^ (( 2P + D C ~ ( 2A;3 + 

(3.3) 

with inverse given by 
1 



■it 



VT 



-.(k(cy-x) -2fcVl-c 2 ) 



,(k{y-cx) -2k 4 VT^), withc = e- i/fc2 



(3.4) 

So, f n (c;x,y) satisfies equations ()2.8j) . which written out, has the form 
'' J ^ {(? - + 2cnx - (c 2 + l)ny 



dx 2 dy \ dxdc 

d3 f f.,2 , a i)2 f o .... , 2 



c - 1) -5-5- + 2cnj/ - (c 2 + l)raa; 



dxdy 2 \ dydc 



2 2 „ 3 2 / A / a 3 / 3 s / 



Using ()3.3|) and definition (2.7), one defines, at a first stage, a new function 

g n (c(t);u,v) := f n (c(t); x(c(t);u,v),y(c(t);u,v)) 

x(c; u, v) — cy(c; u, v) cx(c; u, v) — y(c\ u, v) 



F n c(t); 



c 2 - 1 c 2 - 1 



one expresses the (x, y, c)-partials of / in terms of (u, v, c)-partials of g, using 
the inverse map (|3.4|) . e.g., 

d 3 f ( k \ 3 / <9 3 g <9 3 g 2 d 3 g 2 <9 3 g 



) + (2c 2 + l)-^- + c(c 2 + 2) 



dx 2 dy \ yj l — c 2 / V ^ 3 < ^' y3 du 2 dv dudv 

etc... .One substitutes these expressions in the differential equation ()3.5|) 
for f(c;x,y), yielding a differential equation in g(c;u, v), with coefficients 
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depending on u,v,c and k. After division by k 4 , this differential equation is 
a quadratic polynomial in fc 4 , 



k 8 (c-l) 2 2 



d 3 g cftg 



du 2 dv dudv 2 
d 3 g 



du 2 dv 
d 3 g 



2(c 2 - c + 1)^4 - 4c?4 + v(c 2 + 1) - 2cm 



2(c 2 -c+l 



9m 2 
d 2 g 



dv 2 



+kU 



dudv 2 

d 3 g d 2 g d 3 g d 2 g 



^- 4C + "( C2 + 1 )- 2 - 



-2c 



d 3 g 



&g \ d 2 g \ 



du 3 dv 2 dv 3 du 2 ' \du 2 dv dudv 2 J dudv J 
d 2 g ( d 3 g d 3 g* 



-2 (c 2 -c+ 1 



+c(c 2 - 1] 



dudv \du 3 dv 3 
d 3 g d 3 g 



du 2 dc dv 2 dc 



Then taking into account the fact that 



c(t) 



-t/k 2 



t 

T 2 



o 



d_ 

dc 



k 2 d 
~dt'' 



>+.... 



(3.6) 



(3.7) 



the leading term has order k 4 for large k (noting that k 4 (c — l) 2 = t 2 + 
0(l/k 2 )), At a second stage, defining 



h(t\ u, v) := lim g(e 



-t/k\. 



k— >oo 



and using the expansion ()3.7|) for c(t) and the partial d/dc = —{k 2 /c) d/dt, 
the leading term in the expression ()3.6|) has order k A ; no contribution comes 
from the /c°-term. This leading term (multiplied with —1/2) must therefore 
vanish, leading first to equation (jl.lOj) and then setting x = u — v, y = u + v 
in that equation, to equation (1.6) given in Theorem 1. ■ 
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4 Proof of Theorem 1.2 

Before giving the proof of Theorem 1.2, we remind the reader of the con- 
jecture stated in section 1: for any fixed t > 0, x G M, the conditional 
probability satisfies: 

lim P(A(t) >x + z\ A(0) < -z) = 0. (4.1) 

Under this assumption, we prove the following: 

Lemma 4.1 Considering the series for the probability, for large t, 

P(A(0) < u, A(t) <v) = F 2 {u)F 2 {v) (l + |] ^f^j , (4.2) 

the coefficients fi(u,v) have the property 

lim fi(u,v) = lim fi(v,u) = 0, for fixed uei (4.3) 

and 

lim z + x) — 0, for fixed xeR. (4.4) 

Proof: First observe by (5.1) in section 5 that the Airy kernel becomes di- 
agonal when t — > oo. Then the Airy process decouples at oo, and, using the 
stationarity, one is lead to 

lim P(A(0) < u,A(t) <v) = P(A(0) < u)P(A(0) < v) = F 2 (u)F 2 (v). 



t~*oo 



The next terms follow from the PDE (1.6), although it is more convenient 
here to use the form (1.10) of the equation. Considering the following con- 
ditional probability, 

P(A(t) < v | A(0) < u) 

P(A(0) < u,A(t) < v) 
P(A(0) < u) 



= F 2 (v)(l + Y, 

V i>i 
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and setting 

v = z + x, u = —z, 
we have for all t, since lim^oo F 2 (z + x) = 1, and by (4.1) that 
1 = lim P(A(t) <z + x\ A(0) < -z) 

z— >oo 

]imfi(-z, z + x) 



i>l 



implying that 

lim fi(—z, z + x) = 0, for all % > 1. 

z— >oo 

Similarly, letting u — > oo, we have 



1 = lim P(A(t) < v | A(0) < u) = lim 



*(«>(n-E^) 



lim/i(«,v) 

i>l 

Hence 

lim = (4.5) 

and, considering the same argument for the conditional probability P(A(0) < 
u | A(t) < v), 

lim fi(u,v) = 0, 

ending the proof of Lemma 4.1. ■ 



Proof of Theorem 1.2: Putting the log of the expansion (4.2) 
h(t; u, v) = log P(A(0) < u, A{t) < v) 

= \ og F 2 (u) + \ogF 2 (v)+J2^-^ 

= log F 2 (u) + log F 2 {v) + + + ..., 

(4.6) 

in the equation (jl.l()|) . leads to: 
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(i) a leading term of order t, given by 

Chi = 0, (4.7) 

where 

C—(—- (A 8) 

\du dv J dudv' 

The most general solution to (4.7) is given by 

hi(u,v) = ri(u) + r 3 (v) + r 2 (u + v), 
with arbitrary functions ri,r 2 ,r 3 . Hence, 

P(A(0) < u, A(t) <v) = F 2 (u)F 2 (v) ^1 + M^l + . 

with hi(u, v) = fi(u,v) as in (4.2). Applying Lemma 4.1, 

r\{u) + r 3 (oo) + r 2 (oo) = 0, for all uGl, 

implying 

ri(u) = constant = ri(oo), 

and similarly 

r 3 (-u) = constant = r 3 (cxo). 

Therefore, without loss of generality, we may absorb the constants ri(oo) 
and r 3 (oo) in the definition of r 2 {u + v). Hence, from (4.6), 

fi(u,v) = hi(u,v) = r 2 (u + v) 

using (4.5), 

= lim fi(—z, z + x) = r 2 (x) 

z— +oo 

implying that the hi(u, v)-term in the series (4.6) vanishes. 

(ii) One computes that the term h 2 (u,v) in the expansion (4.6) of h(t;u,v) 
satisfies 
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This is the term of order t°, by putting the series (4.6) in the equation (|1.1(J|) . 
The most general solution to (4.9) is 

h 2 (u, v) = g'(u)g'(v) + r\{u) + r 3 (v) + r 2 (u + v). 

Then 

P(A(0) < u, A(t) < v) = e h{t ' u ' v) 

hj(u,v) 

= F 2 {u)F 2 {v)e^ 2 

- F,(u)F 2 ( V )(l + h -^l + ...). 

In view of the explicit formula for the distribution F 2 and the behavior (1.5) 
of q(a) for a /* oo, we have that 

lim g'(u) = lim (logF 2 (u))' 

u — >oo u^oo 

POO 

= lim / q 2 (a)da = 0. 

u — >oo / 
J u 

Hence 

= lim f 2 (u, v) = lim h 2 (u, v) = ri(oo) + r 3 (v) + r 2 (oo), 



showing ri and similarly r3 are constants. Therefore, by absorbing ri(oo) 
and r 3 (oo) into r 2 {u + v), we have 

f 2 (u,v) = h 2 {u,v) = g'(u)g'(v) + r 2 (u + u). 

Again, by the behavior of q(x) at +oo and — oo, 

/"OO poo 

g'{-z)g'{z + x)= q 2 (a)da q 2 {a)da < c^V 22 / 3 . 



Hence 

= lim f 2 (—z, z + x) = r 2 (x) 

z—*oo 

and so 

f 2 (u,v) = h 2 (u,v) = g'(u)g'(v), 
yielding the 1/t 2 term in the series (4.6). 
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(iii) Next, setting 

h(t- u, v) = log P{A{0) < u, A(t) < v) 

q'(u)q'(v) h 3 (u,v) 
= g(u) + g(v) + y{ >p } + + (4.10) 

in the equation (jl. 10)1 . we find for the t~ l term: 

Ch 3 = 0. 

As in (4.7), its most general solution is given by 

h 3 (u, v) = ri (u) + r 3 (v ) + r 2 (u + v). 
By exponentiation of (4.6), we find 

g'(u)g'(y) 



P(A(0) < u,A(t) <v) = F 2 (u)F 2 (v) 1 



f 2 



ri{u) + r 3 (v) + r 2 (u + v) 
+ p + - 



The precise same arguments lead to h 3 (u,v) = 0. 
(iv) So, at the next stage, we have 



with 



u, u) = + g(v) + 1 g 1 ; + -^f^ + . . . (4.11) 



/ 4 (u, u) = fc 4 (u, «) + 2^( u > v ) = h *( u > v ) + \9\ufg\v) 2 . (4.12) 
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Setting the series (4.11) in the equation (jl.lUj) . we find for the t 2 term: 

£tlA = 



^9_(^9_\ _^9_f^9_\ \ { 9*9 d 3 g ( dg dg 
du 3 \ dv 2 I dv 3 \ du 2 1 j du 3 dv 3 \ du dv 



+ 



1 fd^g_d_ (dg_ 

2 \du7dv [dv 



d 4 g d ( dg^ 



dv 4 du \ du 



f d^g_^g_ + d*g d?g \ _ + ^ ( d^g_dg_ _ d 3 gdg\ 
\du? dv 2 dv 3 du 2 J \du 3 dv dv 3 du J 



oo 

2(2q(u)q'(u)(q(v)q'(v) + l)-q(u)q"(u)q 2 (v)-(q'(u)) 2 q 2 (v)) / q 2 
+ 2q(u)(q(u)q'(v)q"(v) + q'(u)q(v)q"(v) - 2q(u)q 3 (v )q'(v 



— same with u <-» v. (4-13) 

POO 

This latter is an expression in q(u), q{v) and its derivatives and in / q 2 (a)da 

J u 

i>oo 

and / q 2 (a)da. It is obtained by substituting in the previous expression 



g(u) = (u — a)q (a)da 

J u 

and the Painleve II differential equation for q(u), 

u q{u) = q"{u) — 2q(u) 3 , 
in order to eliminate the explicit appearance of u and v. 
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Now introducing 3 

poo 

g(u) = I (u — a)q 2 (a)da 

J u 

POO 

gi(u) = I (u — a)q' 2 (a)da 

J u 

POO 

92(11) = / (u — a)q 4 (a)da : 

J u 

the most general solution to equation (4.13) is given, modulo the null-space 
of C, by 

h,(u, v) = \ {g"(u)g'(vf + g"(v)g'(uf + g"(u)g"(v)) 
+ g'(u)(2g(v) + g' 1 (v)-g' 2 (v)^ 
+ g'(v)(2g(u) + g' 1 (u)-g > 2 (u)^ 



q 2 (u) q \a)da 



V 



+ [ q 2 (a)da f h(v - a)q 2 (a) + q' 2 (a) - q 4 (a))da 

J U J V 

+ same with u <-> v. (4-14) 
This form, together with (4.12), implies for the function / 4 (m, f): 

fi(u,v) = h A (u, v) + ^g'(u) 2 g'(v) 2 + n(u) + r 3 (v) + r 2 (u + v) 
= ^2ai(u)bi(v) + r ± (u) + r 3 (v) + r 2 (u + v). 

i 

Using the asymptotics of q(u), one finds 



3 Notc 



and 



/>oo 

g'(u) = / q 2 (a)da, g"(u) = -q 2 (u), 

J U 

g[(u) = I q' 2 (a)da, g' 2 (u) = I q i (a)da. 

J U J U 
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di(u),bi(u) < ce u u — ► oo, 



< c\u\ 3 u — > — oo, 



and so, by the same argument, 

rx{u) = r 2 (u) = r 3 (u) = 0. 

Therefore, we have 

U(u,v) = h A (u,v) + ^g'(u) 2 g'{v) 2 . 

with /i 4 (m, v) as in (4.14), thus yielding the formula ()1.8|) . 

Finally, to prove formula (jl.9j) . we compute, after integration by parts, 

E(A(Q)A(t)) = // uv— — P(A(0) < u,A(t) < v)du dv 
\ / J J R 2 au dv 

oo poo 

uF2(u)du / vF2(v)dv 



oo J — oo 

oo />oo 



roo /"OO 
' •/— OO J —OO 



+ O 



where 

c:= ^<&(m,v) + <&(?;, = 2 // &(u,v)dudv, 

thus ending the proof of Theorem 1.2. 
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5 The extended Airy kernel 



The joint probabilities for the Airy process can also be expressed in terms 
of the Fredholm determinant of a matrix kernel, the so-called extended Airy 
kernel (considered in [I] and namely 



P{A{t x ) < u 1 ,...,A(t m ) < u m ) 
= det (i - z{k ij ) 1 <i >j < m 



I 2 = 1 



1 + V] ; , / det(k iris (x r ,x s ) ) dx x ...dx k 

/ — ' K\ Jmk V / Kr,s<k 



fc=l l<ii<...<ifc<m 

with 



2 = 1 



Kij{x,y) := X[ui,oo){x) (y) (5.i) 

and 

/ °° e- z( - u ^Ai(x + z)M(y + z)dz, if t t > tj 
- J ^ e- z{u ~^Ai(x + z)k\{y + z)dz, if U < tj , 



Kij(x,y) :-- 



Ai(x) being the Airy function. So the Fredholm determinant above is also a 
solution of the PDE of Theorem 2.1 for m = 2. 



6 Appendix: remarks about the conjecture 

Consider the Dyson Brownian motion (\i(t), . . . , \ n (t)) and the correspond- 
ing Ornstein-Uhlenbeck process on the matrix B. Then, using the change of 
variables 

Bi 

Mi = , 

and further M2 1— > M := M2 — cM\ in the M2-integrals below and noting that 
max(sp Mi) < — z and max(sp M 2 ) > a imply max(sp (M 2 — cM\) > a + cz, 
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we have for the conditional probability, the following inequality: 

P (X n (t) > a | A n (0) < -z) 

J max(sp Mi)<-z J max(sp M 2 )>a 

J max(sp Mi)<-z J M 2 GW„ 

/ dM 1 e-ia-^^M? I dMe -i T rM* 

^ J max(sp M-i_)<-z J max(sp M)>a+cz 

[ dM ie -^-^ M ? [ dM e~^ M2 

J max(sp Mi)<-z J MeH n 

= P(X n (t)>a + cz), 
implying 

lim P (X n (t) > a I A„(0) < -z) = 0, 

and a fortiori, 

lim P (X n (t) > x + z I A n (0) < -z) = 0. 

2— >oo 

It is unclear why this limit remains valid when n — > 00, using the Airy scaling 
(1.3). But the extended Airy kernel (5.1) seems to indicate the conjecture is 
valid. 
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